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Abstract 

In  this  work  the  effect  of  gas  flow-field  design  in  the  bipolar/end  plates  on  the  steady  and  transient  state  performance  of  the  polymer 
electrolyte  membrane  fuel  cell  (PEMFC)  is  presented.  Simulations  were  performed  with  different  flow-field  designs,  viz.  (1)  serpentine;  (2) 
parallel;  (3)  multi-parallel;  and  (4)  discontinuous.  The  steady-state  voltage  at  fixed  current  density  of  5000  A  m~2  was  highest  for  discontinuous 
design.  For  studying  the  transient  response,  the  average  current  density  was  increased  suddenly  from  5000  to  8000  Am-2.  It  was  seen  that 
when  the  load  level  was  increased,  the  voltage  level  suddenly  dropped  and  then  with  time  leveled  off  to  a  value  slightly  higher  than  the  dropped 
value.  This  time  for  serpentine,  parallel,  multi-parallel  and  discontinuous  flow-fields  were  9.5,  7.5,  8.0  and  16.5  s,  respectively.  While  it  was 
seen  that  the  steady- state  performance  of  the  discontinuous  type  of  design  was  the  maximum,  its  transient  response  was  slow.  On  the  other 
hand  in  case  of  parallel  type  of  design  the  steady-state  performance  was  low,  but  the  transient  response  was  high.  The  multi-parallel  design 
offers  a  unique  advantage  of  both  of  these  properties,  viz.  high  steady-state  performance  with  good  transient  response,  and  therefore  should 
perform  better  than  the  other  designs  chosen  in  this  study. 

©  2005  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

The  polymer  electrolyte  membrane  fuel  cell  (PEMFC)  is 
seen  as  a  system  of  choice  for  automotive  applications  due 
to  several  advantages  that  this  type  of  fuel  cell  offers  over 
the  other  types  (alkaline  fuel  cell,  phosphoric  acid  fuel  cell, 
molten  carbonate  fuel  cell  and  solid  oxide  fuel  cell).  PEM- 
FCs  operate  at  lower  temperature  (~80°C),  thereby  allow¬ 
ing  for  faster  startups  and  immediate  response  to  changes 
in  demand  of  power.  Also  these  systems  are  environmen¬ 
tally  friendly  and  hence  suited  best  for  automobile  applica¬ 
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tions.The  PEMFC  consists  of  an  anode,  cathode  and  a  proton 
exchange  membrane  (electrolyte).  The  electrolyte  membrane 
is  sandwiched  in  between  anode  and  cathode.  Pure  hydrogen 
or  hydrogen  rich  reformate  gas  is  passed  onto  the  anode  side, 
wherein  the  hydrogen  breaks  into  protons  (H+)  and  electrons 
(e-).  These  H+  ions  pass  through  the  electrolyte  (protonic 
conductor)  onto  the  cathode  side,  while  the  electron  follows 
through  the  outer  circuit  generating  electricity.  Air  or  oxygen 
is  the  oxidant  on  the  cathode  side.  The  oxygen  molecule  com¬ 
bines  with  these  H+  ions  and  electrons  on  the  cathode  side  to 
generate  the  by-product  water  and  heat.  The  electrochemical 
reactions  for  the  process  can  be  represented  as: 

anodic:  H2  — >  2H++2e_  (1) 
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Nomenclature 

ak 

activity  of  water  at  interface  k 

Acv 

specific  area  of  control  volume  (m2  m-3) 

D,j 

binary  mass  diffusion  coefficient  of  compo¬ 
nent  i  in  component  j  (m2  s-1) 

E 

theoretical  cell  voltage  (V) 

F 

Faradays  constant  (C  eq-1) 

^stk 

compaction  force  of  the  stack  (N  m-2) 

Iq 

exchange  current  density  for  the  oxygen  reac¬ 
tion  (100  Am-2) 

I 

current  density  (A  m-2) 

/l 

limiting  current  density  (Am-2) 

h 

diffusion  flux  of  species  k  (kg  m-2  s-1) 

kp 

permeability  of  water  in  membrane  (m2) 

mk 

mass  fraction  of  species  k 

Mk 

molecular  weight  of  species  k  (kg) 

Mm  Ary 

equivalent  weight  of  dry  membrane 
(kg  mol-1) 

Pk 

partial  pressure  of  species  k  (Pa) 

P 

pressure  (Pa) 

psat 

^w,k 

Saturation  vapor  pressure  of  water  at  ‘k’  inter¬ 
face 

R 

universal  gas  constant  (J  mol-1  K-1) 

^CR 

contact  resistance  of  material  k  (ohms) 

^mcmb 

resistance  of  membrane  (ohms) 

sk 

source  term  for  species  k  (kgm-3  s-1) 

hn 

polymer  membrane  thickness  (m) 

T 

temperature  (K) 

V 

velocity  vector  (ms-1) 

V 

actual  cell  voltage  (V) 

Greek  symbols 

a 

net  water  flux  per  H+ 

p 

permeability  of  the  electrodes  (m2) 

7?  act 

activation  polarization  (V) 

0  ohm 

ohmic  polarization  (V) 

0  cone 

concentration  polarization  (V) 

P 

density  (kgm-3) 

Pm,  dry 

density  of  dry  electrolyte  membrane  (kgm-3) 

membrane  conductivity  (mhom-1) 

Pk 

viscosity  of  species  k  (Pas-1) 

cathodic:  l/202  +  2H+  +  2e  — >  H2O  (2) 


net-reaction  :  H2  +  1  /2O2 

— >  H2O  +  Heat  +  Electricity  (3) 

In  spite  of  the  presence  of  fuel  cell  technology  for  over  150 
years,  some  technical  issues  exist  that  prohibit  the  wide  scale 
commercialization  of  these  fuel  cells.  Some  of  these  issues 
related  in  particular  to  PEMFC  are:  (1)  high  cost  of  PEMFC 
components;  (2)  catalyst  poisoning  by  CO  concentration  in 
the  fuel;  (3)  water  management  in  large  stacks;  (4)  proton 


exchange  membrane  operating  temperature  limits;  and  (5) 
cell  stack  life.  The  biggest  challenge  of  these  is  the  reduction 
in  cost  of  the  fuel  cell  stack,  which  is  mainly  ascribed  to  the 
bipolar/end  plates  and  precious  metal  catalysts.  Additionally, 
the  cost  can  be  decreased  by  the  optimization  of  the  fuel  cell 
process.  Several  models  have  been  developed  for  optimiza¬ 
tion  of  the  electrochemical  kinetics  of  the  process.  However, 
little  attention  has  been  on  the  optimization  of  the  flow-field 
design  in  the  bipolar/end  plates  of  the  fuel  cell  [1-3].  For 
example,  Kumar  and  Reddy  showed  in  their  work  that  opti¬ 
mum  dimensions  of  the  channels  in  the  flow-field  exist  where 
the  fuel  cell  gives  the  maximum  performance  [3] .  Most  of  the 
existing  models  are  developed  under  steady-state  conditions 
and  assume  the  flow-field  to  be  serpentine  type.  This  work 
therefore  concentrates  on  the  effect  of  the  flow-field  design 
on  the  steady  state  and  the  transient  behavior  of  the  fuel  cell. 
Several  flow-field  designs,  viz.  (1)  serpentine;  (2)  parallel;  (3) 
multi-parallel;  and  (4)  discontinuous  were  studied  for  their 
steady  state  and  transient  performance  of  the  fuel  cell. 

2.  Model  development 

The  current  model  is  a  three-dimensional  single-phase 
transient  isothermal  numerical  mass-transfer  unified  model 
for  PEMFC.  The  steady-state  performance  is  just  a  special 
case  where  all  the  time  dependent  terms  in  the  governing 
transport  equations  would  be  nullified.  Numerical  predictions 
of  the  voltage  would  be  made  for  changing  load  (current) 
levels.  The  transient  behavior  would  be  studied  for  differ¬ 
ent  flow-field  designs  (serpentine,  parallel,  multi-parallel  and 
discontinuous)  in  bipolar/end  plates  of  the  fuel  cell.  Also,  the 
steady-state  performance  would  be  calculated  for  each  of  the 
designs  and  results  would  be  compared.  The  overall  perfor¬ 
mance  of  the  flow-field  design  would  be  judged  based  on 
the  steady  state  and  transient  behavior  of  the  fuel  cell.  This 
would  help  in  finding  the  best  possible  flow-field  design  for 
the  bipolar/end  plates  of  the  fuel  cell.  The  detailed  step-wise 
development  of  the  unified  model  is  discussed  in  following 
sections. 

2.7.  Problem  domain 

The  simulation  domain  consists  of  cathode  and  anode 
flow-fields  corresponding  to  the  cathode  and  anode  sides  of 
the  fuel  cell,  respectively,  separated  by  the  membrane  elec¬ 
trode  assembly  (MEA)  in  between.  The  MEA  consists  of 
cathode  and  anode  porous  gas  diffusion  layer  and  catalyst 
layer,  on  each  side  of  the  membrane.  The  active  area  of  the 
domain  was  a  square  of  side  7.62  cm  (3  in.).  Fig.  1  shows  the 
detailed  cross-section  ( X-Z  section)  of  the  problem  domain 
for  this  model.  Each  of  the  layers  is  marked  by  a  separate 
zone.  These  are:  cathode  flow-field;  cathode  diffusion  layer; 
cathode  catalyst  layer  on  the  cathode  side;  anode  flow-field; 
anode  diffusion  layer;  and  anode  catalyst  layer  on  the  anode 
side  of  the  fuel  cell.  The  cathode  and  anode  flow-field  were 


266 


A.  Kumar,  R.G.  Reddy  /  Journal  of  Power  Sources  155  (2006)  264-271 


+Z  axis 

n 


z„+ 


z  =z  - 

*"09  10 


7  =  7  - 

*"07  *"08 

7  =7  - 

*"05  *"06 

7  =  Z  - 

*"03  *"04 


z  =z  - 
*"0 1  *"02 


z  -- 


"00 


Anode  Flow-Field  Channel  (1/16”) 

Anode  Diffusion  Layer  (1/64”) 

Anode  Catalyst  Layer  ( 1/320”) 
Electrolyte  Membrane 
Cathode  Catalyst  Layer  (1/320”) 

Cathode  Diffusion  Layer  (1/64”) 

Cathode  Flow-Field  Channel  ( 1/16”) 


Diagram  not  to  scale 
•  AH  Dimension  in  Inches 


X 

— ► 


Y 


Fig.  1.  Schematic  of  polymer  electrolyte  membrane  fuel  cell  cross-section  (X-Z  section)  showing  different  zones  and  transport  of  species  across  these  zones. 


essentially  the  machined  channels  in  the  bipolar/end  plates. 
These  help  to  distribute  the  reaction  gases  uniformly  over 
the  reaction  surface.  In  this  work  four  different  flow-field 
designs  were  studied.  Fig.  2(a-d)  shows  the  schematic  (X-Y 
section)  for  all  of  these  flow-field  designs.  The  dimension 
of  the  channel  width,  land  width,  and  channel  depth  was 
1.587  mm  (1/16  in.).  The  cathode  and  anode  gas  diffusion 
layers  were  the  electrodes  made  with  porous  carbon  cloth 
that  helps  to  diffuse  the  reactant  gases  from  the  bipolar/end 
plate  flow-field  channels  towards  the  reaction  catalyst  layer, 
and  also  diffusion  of  byproduct  water  from  the  reaction  site 


(C) 


back  to  the  flow-field  channels.  The  catalyst  layer  on  both  the 
cathode  and  anode  sides  consists  of  platinum  coating,  and  is 
essentially  where  all  the  electrochemical  reactions  take  place. 

We  consider  four  different  species  in  our  model,  viz. 
hydrogen  (H2),  oxygen  (O2),  nitrogen  (N2),  and  water  (vapor) 
(H2O).  Each  of  the  species  flux  was  tracked  in  the  domain 
through  the  control  volume  technique,  which  was  used  in 
solving  the  model.  For  example,  consider  the  hydrogen 
species.  The  hydrogen  from  the  anode  flow-field  channels 
is  transported  through  the  anode  diffusion  layer  towards 
the  anode  catalyst  layer.  The  catalyst  action  dissociates  the 


Fig.  2.  Schematic  of  (a)  serpentine;  (b)  parallel;  (c)  multi-parallel;  and  (d)  discontinuous  flow-field  designs  used  in  simulations. 
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hydrogen  molecule  to  protons  and  electrons  according  to  the 
Eq.  (1).  The  protons  generated  above  will  transfer  through  the 
electrolyte  membrane  onto  the  cathode  side,  wherein  they 
combine  with  oxygen  to  generate  water.  So  there  is  a  con¬ 
sumption  of  hydrogen  and  oxygen  species  on  the  anode  and 
the  cathode  side,  respectively,  and  generation  of  water  species 
on  the  cathode  side.  Similarly  other  species  were  tracked  in 
the  domain. 

2.2.  Model  assumptions 

The  following  assumptions  were  used  in  developing  the 
unified  model  for  PEMFC. 

(1)  Stationary  conditions  exist  in  the  single-cell  PEMFC 
stack.  Also,  the  effect  of  gravity  was  neglected. 

(2)  Based  on  the  Reynolds  number  calculation  (Re  ~  500) 
the  flow  in  the  fuel  cell  domain  is  laminar.  Hence,  all 
the  transport  equations  were  formulated  for  the  laminar 
regime  behavior. 

(3)  The  proton  exchange  membrane  was  considered  as  a 
solid  wall  in  the  development  of  the  model,  and  trans¬ 
port  of  species  (protons  and  water  molecules)  across  the 
membrane  was  taken  into  account  through  the  incorpo¬ 
ration  of  source  terms  in  the  transport  equations. 

(4)  The  gas  diffusion  layer  and  the  catalyst  layer  on  both 
the  cathode  and  the  anode  sides  were  assumed  to  be 
isotropic  with  a  permeability  of  10“ 12  m2. 

(5)  Four  reactant  species,  viz.  H2, 02,  N2,  and  H2O  (vapor) 
were  considered. 

(6)  Bipolar/end  plates  fabricated  out  of  SS-316  alloy  were 
used  in  the  development  of  the  model. 

(7)  The  volume  of  the  liquid  H2O  was  assumed  to  be 
negligible  in  the  domain.  Due  to  the  high  thermal  con¬ 
ductivity  of  metallic  bipolar/end  plates,  it  was  assumed 
that  whatever  liquid  water  droplets  are  present  would 
be  only  on  the  walls  of  the  channels  in  the  flow-field 
of  these  plates.  Therefore,  it  can  be  assumed  that  there 
is  no  restriction  to  the  flow  of  gases  in  the  thin  chan¬ 
nels  of  the  flow-field  due  to  the  presence  of  liquid  water 
species. 

(8)  Due  to  pressure  variation  and  electrochemical  reactions, 
the  density  of  the  gaseous  mixture  varies  from  loca¬ 
tion  to  location  in  the  domain.  Hence,  compressible  gas 
technique  (ideal  gas  behavior)  was  used  to  determine 
the  density  of  the  gas  mixture  at  every  control  volume 
element  in  the  domain. 

(9)  Gas  mixture  viscosity  was  calculated  using  mass- 
weighted-mixing-law. 

(10)  The  species  binary  diffusion  coefficient  was  calculated 
using  kinetic  theory  of  gases  [4]. 

(11)  Since  we  are  using  metallic  bipolar/end  plates  (without 
surface  modification),  there  is  an  appreciable  level  of 
contact  resistance  from  these  plates.  To  take  these  losses 
into  account,  we  assume  that  the  oxide  film  thickness 
on  these  plates  to  be  constant,  and  that  the  contact  resis¬ 


tance  of  the  bipolar/end  plates  with  the  electrodes  was 
a  function  of  compaction  pressure  of  stack. 

2.3.  Governing  equations 

The  governing  equations  for  this  unified  model  can  be 
divided  into  three  different  categories,  viz.  (1)  transport  equa¬ 
tions;  (2)  thermodynamic  equations;  and  (3)  membrane  prop¬ 
erties  equations. 


2.3.1.  Transport  equations 

These  include  equations  of  conservation  of  mass,  momen¬ 
tum  and  species  transport  in  the  domain.  These  equations  are 
implemented  zone-wise  in  the  model  domain.  The  equation 
for  conservation  of  mass  in  general  form  can  be  described 
as: 

dp 

i-  +  V  •  (pv)  =  (4) 

where  p  is  the  density  of  the  gaseous  mixture;  t  the  time;  v  the 
velocity  vector;  and  Sm  is  the  source  term  for  the  continuous 
phase.  This  source  term  corresponds  to  the  production  or 
consumption  of  reaction  species  in  a  particular  zone.  Since  all 
the  electrochemical  reactions  take  place  only  in  the  cathode 
and  anode  catalyst  layers,  the  value  of  Sm  is  only  non-zero 
for  these  layers  and  is  given  by: 


5m,  c  —  5o2  T  3w  c 
5m, a  =  5h2  T  5w?a 


(5) 

(6) 


where  the  subscripts  ‘c’  and  ‘a’  refer  to  the  cathode  and  anode 
sides,  respectively,  and  Sx  (x  =  H2,  O2,  w  (water))  represent 
the  source  terms  for  species  x  in  the  catalyst  layer.  The  density, 
p  mixture  is  calculated  using  the  compressible  gas  technique 
(ideal-gas  behavior)  and  is  given  for  multi-component  system 
as  [4]: 


P  = 


Pop  P 

m; 
Mi 


RTJ2 


where  PGp  is  the  operating  pressure;  P  the  local  relative 
(or  gauge)  pressure,  m;  the  mass-fraction  of  species  /,  and 
Mi  is  the  molecular  weight  of  species  i. The  conservation 
of  momentum  equation  in  general  form  can  be  described 
as: 


— (pv)  +  v.(pw) 

Ot 


— V P  +  V  •  (r)  +  Sp  +  5m,k^ 


where  r  is  the  stress  tensor  and  Sv  is  the  momentum  source 
term.  For  laminar  flow  through  porous  media,  the  momentum 
source  term  Sv  is  proportional  to  velocity  and  is  given  by 
Darcy’s  law: 
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where  fi  is  the  viscosity  of  the  gaseous  mixture  given  by 
mass- weighted  mixing  law: 

ii = t  miii,  (io) 

i 


and  fi  is  the  permeability  of  the  medium.  The  value  of  /3  was 
taken  to  be  10“ 12  m2  for  the  electrodes  and  catalyst  layer 
The  species  transport  equations  were  written  in  general 
form  for  each  of  the  species  H2,  O2,  N?,  and  H2O  vapor.  On 
the  anode  side,  the  species  transport  equations  are: 

d 

—(pmu2)  +  V(pnmH2)  =  -V  •  7h2  +  5r2  (H) 

ot 

d 

—  (p/fZw,a)  T  V(piWZw?a)  —  —V  •  Jw  a  +  ^w,a  (12) 

ot 


and  for  the  cathode  side  the  species  transport  equations  of 
O2,  N2,  and  H2O  are: 

d 

—  (pmo2)  +  V(pnmo2)  =  -V  •  J0l  +  So 2  (13) 

ot 

d 

—  ( pmNl )  +  V(pvmN2)  =  -  V  •  JNl  (14) 

ot 

d 

—  (pmwc)  T  V (pvmwc)  —  —V  •  /w?c  T  Sw,c  (15) 

ot 

where  Ji  is  the  diffusion  flux  of  species  i.  The  diffusion  flux 
of  species  i  is  given  by  Maxwell  relationship  as  [4] : 


N- 1 

7  =  -  f>D'-  iVm  i  (16) 

7=1 


where  Dq  is  the  binary  diffusion  coefficient  of  species  i  in 
species  j,  and  N  is  the  total  number  of  species  in  the  mixture. 
The  binary  diffusion  coefficients  (Pq)  were  calculated  using 
the  kinetic  theory  of  gases  [4] .  The  source  terms  for  each  of 
the  species  transport  equations  exist  only  in  their  respective 
catalyst  layer  and  are  given  by: 


’w,a 


I(x,  y)~ 

IF 

a(x,  y) 

F 

I(x,  }’)" 

4  F 

Mr,  A 


2^cv 


(17) 


I(x,  y)Mn2oA 

CV 


(18) 


Mq,  A 


2^cv 


(19) 


1  +  2a(x,  y) 
2  F 


I(x,  y)Mu2oA  CV 


(20) 


where  I(x ,  y)  is  the  local  current  density;  F  the  Faradays  con¬ 
stant;  a(x ,  y)  the  local  net  water  transfer  coefficient  per  proton; 
and  Acv  is  the  specific  surface  area  of  control  volume  element 
in  the  respective  zone  of  the  domain. 


2.3.2.  Thermodynamic  equations 

The  Nemst  equation  for  the  H2/O2  fuel  cell,  using  litera¬ 
ture  values  for  the  standard- state  entropy  can  be  written  as: 

E  =  1.229  -  (8.5  x  10“4)(T  -  298.15) 

+(4.308  x  10-5)T(ln(/4)  +  \  In (/>*2))  (21) 


where  E  (V)  is  the  ideal  equilibrium  potential,  T  (K)  is  the  cell 
temperature  and  Pf  (atm)  is  the  partial  pressures  of  species 
i  at  the  reaction  interface.  However,  the  actual  potential  of 
the  cell  decreases  when  the  load  is  drawn  from  the  cell.  This 
is  due  to  irreversible  losses  in  the  cell  which  include  viz.  (1) 
activation  polarization,  ^act;  (2)  ohmic  polarization,  ?7ohm; 
and  (3)  concentration  polarization,  qconc.  These  losses  are 
represented  by  standard  expressions  [5]  and  are  given  in  Eqs. 
(22)— (24). 


_  RT  ( I(x,  y) 

7?  act  —  hi 


F 


Io 


+>hm  =  Kx.  y)Rcen(x,  y) 

_RT  f  I(x,  y) 
Oconc  —  In  (  1 


nF 


Il(x ,  y) 


(22) 


(23) 


(24) 


where  Io  is  the  exchange  current  density;  RCQ\\  the  ohmic 
resistance  offered  by  the  cell  stack;  and  7l  the  limiting  current 
density.  The  actual  cell  voltage  (V)  is  then,  theoretical  cell 
voltage  ( E )  minus  the  polarization  losses. 


V  —  E  Oac  t  77ohm  Oconc 


(25) 


The  ohmic  resistance,  Rce n  of  the  cell  stack  is  a  combination 
of  the  resistance  to  flow  of  electrons  in  cell  components  (typ¬ 
ically  bipolar  plates)  and  the  resistance  to  flow  of  ions  (H+) 
in  the  electrolyte  and  is  given  by: 

^ceiiC*,  y)  =  Rcr(x,  y )  +  ^membC*,  y)  (26) 


where  Rqr(x,  y)  is  the  contact  resistance  of  the  stack  and 
^memb  is  the  resistance  offered  by  the  electrolyte  membrane 
to  the  flow  of  H+  ions.  The  contact  resistance  depends  upon 
the  type  of  material  for  the  bipolar/end  plates  and  compaction 
force  of  the  stack.  We  have  shown  in  our  previous  work  that 
metallic  bipolar/end  plates  have  several  advantages  over  the 
conventionally  used  graphite  plates  [6].  Hence,  we  decided 
to  develop  this  model  using  SS-3 16  bipolar/end  plates.  How¬ 
ever,  the  contact  resistances  while  using  these  materials  (with 
no  surface  modification)  can  be  significant  due  to  the  forma¬ 
tion  of  thin  oxide  layer  on  the  surface.  Consequently,  these 
contact  losses  should  be  considered  in  the  modeling  work. 
A  relationship  between  compaction  force  and  contact  resis¬ 
tance  for  SS-3 16  was  developed  by  curve  fitting  the  data 
from  Hodgson  et  al.  [7].  The  equation  for  contact  resistance 
for  SS-3 16  is  given  by: 


r>SS316  _  0  ^  77-0.9231 
^CR  —  Z.OZTstk 


(27) 


where  ^^R316  contact  resistance  (fi  m2)  and  Fs±  is  the 
compaction  pressure  of  the  stack  in  N  m-2.  The  resistance  to 
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Table  1 

Equations  for  local  membrane  properties  [8-12] 

Activity  of  water  (a)  on  the 
membrane-electrode  interface 

Saturation  vapor  pressure  of  water  at  ‘k’ 
interface 

Water  content  in  membrane  (k  =  a) 

Electro-osmotic  drag  coefficient 

Diffusion  coefficient  of  water  in  the 
membrane  at  temperature  Tc  (K) 


Water  concentration  at  membrane-electrode 
interface 


Membrane  conductivity 


Net  water  transfer  coefficient  in  membrane 
neglecting  the  convection  term 


«k  = 


*w,k  P 
nsat 
w,k 


togio^k)  =  2-95  x  10“z7k  -9.18  X  ur*T4  +  1.44  x  W~‘T>  -2.18 


5  t^2 
k 


-7  ^3 
k 


A  =  0.043  +  17.81a);  —  39.85ak  +  36 ak,  0  <  ak  <  1;  A  =  14  +  1.4 (ak  —  1), 

1  <  «k  <  3 

tii  =  0.0029A2  +  0.05  A  —  3.4  x  10“ 19 

£>w  =  A  exp  (2416  x  (4-  -  ^))  ,  7b  =  303K;  D,  =  10-10,  X  <  2,  D,=10-10, 

A  <2;  Dj  =10“ 10,  (1  +2(A  —  2)),  2<A<3;  D,  =  10  10(3  —  1.67(A  —  3)),  3<\<4.5; 
Dj  =  1.25x  10-10,  A  >4.5 


CWik  =  f”,dly  (0.043  +  17.8ak  -  39.85ak  +  36a3k)  ak  <  1 


^m.dry 

C„,k=  f^(14  +  1.4(ak- 1)) 


j  Pm, dry  —  2000  kg  m  ^ , 


AJ  m.dry 

Mm,dry  =  l.lkgmor1 

^m.dry 


>  1 


crm  =  (0.514^^Cw,a-  0.326]  exp  [l268  (T  -  ^)]  ;  T0  =  303K,  resurface 


Pm,  dry 

temperature  at  anode 


a  —  nx\ 


i  ^  (Cw,c  Cw,a) 
7 


flow  of  ions  (H+)  in  the  electrolyte  membrane  depends  on  the 
proton  conductivity  of  the  membrane  (am)  and  the  thickness 
of  the  membrane  (fm),  and  is  given  by  the  expression: 

tfmemb  =  ^  ,  (28) 


2.3.3.  Membrane  properties  equations 

The  water  content  and  the  flux  of  water  in  the  membrane 
are  important  for  current  density  prediction  in  the  cell  [8,9]. 
The  equations  for  the  electrolyte  membrane  properties  are 
empirical  equations  and  are  obtained  by  experimental  data 
[8-12].  Table  1  gives  a  compiled  list  of  all  the  membrane 
properties  equations  that  were  used  in  the  development  of 
this  model  and  the  detailed  explanation  of  these  equations  are 
given  in  several  references  [8-12].  It  may  be  noted  that  all  the 
properties  referred  to  in  this  section  are  the  local  properties 
values  on  the  X-Y  reaction  interface,  and  the  suffix  (v,  y)  has 
been  taken  out  just  for  mere  simplicity. 


2.4.  Solution  strategy 

The  geometry  of  the  domain  was  defined  in  Gambit®  soft¬ 
ware.  The  geometry  consists  of  cathode  and  anode  flow-fields 
and  membrane  electrode  assembly  (ME A)  sandwiched  in 
between.  For  meshing  the  domain,  hexahedron  volume  ele¬ 
ment  with  20  nodes  per  element  and  map  scheme  was  chosen. 
Total  number  of  the  volume  elements  with  the  chosen  mesh 
size  was  around  0.4  million.  The  geometry  was  then  exported 
to  Fluent®  6.0  software.  The  formulated  transport  equations 
were  solved  using  the  simple  algorithm  [4]  for  each  of  the 
control  volume  element  in  the  domain.  Second  order  implicit 
method  was  used  for  unsteady- state  formulation.  Thermody¬ 
namic  and  membrane  properties  equations  were  incorporated 
through  the  use  of  user  defined  function  (UDFs)  codes  writ¬ 
ten  in  C  language.  All  the  equations  were  solved  and  the 


parameters  (membrane  properties,  local  current  density,  etc.) 
were  updated  at  the  start  of  every  iteration  using  the  adjust 
function.  A  grid  adaption  technique  was  used  to  obtain  a 
grid  independent  solution.  Because  of  the  memory  and  time 
requirements  for  the  iteration  process,  the  Fluent®  and  the 
Gambit®  softwares  were  run  on  Dell®  Precision®  worksta¬ 
tion  (Intel®  Xeon®  2.2  GHz,  1.5  GB  SDRAM). 


3.  Results  and  discussions 

Simulations  were  done  for  different  flow-field  designs  in 
the  bipolar/end  plates  of  the  polymer  electrolyte  membrane 
fuel  cell.  Temperature  and  pressure  for  all  the  simulations 
were  kept  constant  at  350  K  and  202  kPa,  respectively.  Reac¬ 
tant  gases  were  externally  humidified.  On  the  cathode  side, 
humidified  air  consisting  of  21  mass%  of  O2,  70.5  mass%  of 
N2  and  8.5  mass%  of  H2O  vapor  was  used.  The  anode  side 
reactant  gas  consists  of  40  mass%  of  H2  and  60  mass%  of 
H2O  vapor.  The  velocity  inlets  for  the  cathode  and  anode 
sides  were  5.0  and  2.5  m  s-1,  respectively. 

The  steady- state  results  were  first  obtained  by  de¬ 
selecting  the  transient  option  from  the  Fluent®  software. 
De-selecting,  essentially  removes  all  the  time  dependent 
terms  from  the  transport  equations.  The  current  density  for 
these  simulations  was  fixed  at  5000 Am-2.  As  discussed 
earlier,  voltage  level  was  used  to  judge  the  performance  of 
the  fuel  cell.  The  simulated  voltages  for  serpentine,  parallel, 
multi-parallel  and  discontinuous  designs  were  0.66,  0.64, 
0.68  and  0.71  V,  respectively.  These  are  marked  by  dotted 
lines  in  the  steady-state  region  of  Fig.  3.  From  this,  we  can 
conclude  that  the  discontinuous  type  of  flow-field  design 
will  perform  better  than  the  other  three  designs.  This  is 
because  the  discontinuity  of  the  channels  forces  the  gas  into 
the  diffusion  layer,  thereby  making  the  transfer  of  reactant 
gases  in  the  gas  diffusion  layer  from  diffusion  to  diffusion 
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Fig.  3.  Steady  state  and  transient  performance  for  one  cell  PEM  fuel  cell 
stack  using  different  flow-field  designs  in  the  bipolar/end  plates. 

plus  convection  type.  This  increases  the  effective  pressure  of 
the  reactions  species  at  the  reaction  interface.  These  findings 
are  consistent  with  those  of  Watkins  et  al.  [1],  who  proposed 
the  design  as  a  solution  to  the  problem  of  increased  gas 
diffusion. 

The  transient  response  of  the  PEMFCS  was  studied  by 
changing  the  load  level  (current)  from  5000  to  8000  Am-2 
instantaneously.  The  sudden  load  change  was  made  at  time, 
t-  0  s.  It  was  assumed  that  the  fuel  cell  has  a  steady  state 
at  all  times,  f<0s.  The  simulated  voltage  for  the  transient 
region  is  shown  in  Fig.  3.  Initially  when  the  load  level  is 
increased,  the  voltage  suddenly  drops,  and  then  leveled  off 
to  a  value  slightly  higher  than  the  dropped  value.  The  time 
to  level  off  was  measured  from  the  time  when  the  load  level 
was  suddenly  increased  (t  =  0  s)  to  the  time  when  the  change 
in  voltage  is  less  than  0.1  mV.  The  time  for  the  transient 
response  for  serpentine,  parallel,  multi-parallel  and  discon¬ 
tinuous  type  of  flow-designs  was  9.5,  7.5,  8.0,  and  16.5  s, 
respectively.  It  was  seen  that  the  parallel  type  of  design  gave 
the  best  performance,  while  the  transient  performance  from 
the  discontinuous  type  of  design  was  low. 

As  discussed  above,  the  discontinuous  type  of  design 
gives  the  maximum  steady-state  performance,  but  lowest 
transient  performance.  On  the  other  hand  the  parallel  flow- 
field  design  gave  the  lowest  steady- state  performance  but 
highest  transient  performance.  The  multi-parallel  type  of 
flow-field  design  on  the  other  hand  gives  reasonable  good 
values  for  both  the  steady  state  and  transient  performance. 
So,  with  such  flow-field  designs  in  the  bipolar/end  plates, 
the  overall  performance  characteristics  of  the  fuel  cell  can  be 
increased. 

Also  it  is  worth  mentioning  that  the  use  of  design  also 
depends  upon  the  type  of  applications.  If  the  fuel  cell  is 
used  for  stationary  applications  where  steady- state  response 
is  more  important  than  the  transient  response,  discontinuous 
type  of  flow-field  design  will  serve  as  a  better  choice.  How¬ 
ever,  if  the  fuel  cell  were  used  for  automobile  applications, 
where  transient  response  is  important  along  with  the  steady- 


state  performance,  it  would  be  wise  to  use  multi-parallel  type 
of  flow-field  designs. 

Also  it  may  by  noted  that  the  voltage  transient  response 
values  are  in  order  of  0. 1  V,  which  may  appear  to  be  insignifi¬ 
cant.  These  values  should  be  used  as  a  caution,  as  this  depicts 
the  voltage  variation  in  one  cell  only.  In  commercial  cell 
stacks  (~100  cell  or  more),  the  voltage  variations  from  all 
the  cells  would  add  up  resulting  in  an  appreciable  level  of 
transient  voltage  drop.  The  present  model  in  our  further  stud¬ 
ies  would  be  extended  to  study  the  transient  characteristics 
for  the  large  fuel  cell  stacks.  The  use  of  such  models  will  help 
in  the  optimization  of  the  design  process  for  the  flow-field  in 
the  bipolar/end  plates  and  will  give  the  fuel  cell  developers  a 
scope  of  further  improvement  of  the  fuel  cell  technology. 


4.  Conclusions 

A  unified,  three-dimensional,  transient  numerical  mass- 
transfer  single  cell  model  for  polymer  electrolyte  membrane 
fuel  cell  was  developed.  Simulations  were  performed  with 
different  flow-designs  in  the  bipolar  plate.  The  relative  per¬ 
formance  of  the  different  designs  was  judged  from  the  voltage 
levels  at  a  constant  current  density.  The  simulated  steady- state 
voltage  for  discontinuous  type  of  design  at  an  average  cur¬ 
rent  density  of  5000  A  m-2  was  0.7 1  V,  which  was  maximum 
amongst  the  designs  chosen  for  simulations  in  this  study. 
The  transient  performance  of  the  cell  was  studied  when  the 
load  levels  were  varied  from  5000  to  8000 Am-2.  Results 
showed  that  the  transient  response  for  the  parallel  type  of 
design  was  highest,  and  it  took  only  7.5  s  for  the  cell  to  get 
back  to  steady  state  again.  However,  the  steady- state  voltage 
for  parallel  flow-field  is  comparatively  lower  than  the  other 
three  designs.  An  optimized  approach  is  therefore  to  choose 
multi-parallel  design,  where  reasonable  values  of  steady- state 
voltage  at  5000 Am-2  and  transient  response  are  obtained. 
This  model  can  be  extended  to  calculate  the  steady  state  and 
transient  performances  of  the  large  commercial  stacks,  and 
can  help  in  to  improve  the  performance  characteristics  of  the 
fuel  cell  stack. 
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